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Abstract. We consider Frenkel-Kontorova models corresponding to 1 
dimensional quasicrystals. 

We present a KAM theory for quasi-periodic equilibria. The theorem 
presented has an a-posteriori format. We show that, given an approxi- 
mate solution of the equilibrium equation, which satisfies some appro- 
priate non-degeneracy conditions, then, there is a true solution nearby. 
This solution is locally unique. 

Such a-posteriori theorems can be used to validate numerical compu- 
tations and also lead immediately to several consequences a) Existence 
to all orders of perturbative expansion and their convergence b) Boot- 
strap for regularity c) An efficient method to compute the breakdown of 
analyticity. 

Since the system does not admit an easy dynamical formulation, the 
method of proof is based on developing several identities. These identi- 
ties also lead to very efficient algorithms. 

Quasi-periodic solutions, quasicrystals, hull functions, KAM theory 
[2000] 70K43, 52C23, 37A60, 37J40, 82B20 



1. Introduction 

We consider Frenkel-Kontorova models with a quasi-periodic potential. 
In these models, the state of a system is given by a sequence {<?i},gz of num- 
bers and the physical states are selected to be the critical points of a formal 
energy 

(1) =^({M}iez) = ~ ""+^'*^ ~ ^^""-^ 

fieZ 

The critical points of are obtained by taking formal derivatives of 
and setting them to zero ^^({M},ez) = 0. That is 

(2) M„+l + Un-\ - 2Un + V'{Un) = 0. 

In the classical Frenkel-Kontorova model IIFK391 IALD83I1 . V is a pe- 
riodic function V{Un + 1) = y(M„)-we choose units so that the period is 
normalized to 1 . 
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In this paper, however, we choose V to be a quasi-periodic function 

V{e) = ViOa) 

where V : T"^ ^ R. and 6 e R, a e Mf' will be an irrational vector, that is 
k - a when k eZ'' - {0} where d >2. Later, we will also assume further 
non-resonance conditions. One example to keep in mind could be 

V(e) = a ■ smilnO) + b ■ sm(2n ^flO) 

where a, 6 R. 

Our goal in this paper is to prove an analogue of KAM theorem. See 
Theorem [m [T^l Under some appropriate conditions, we show that there 
are smooth families of quasi-periodic solutions. See Section [TTT] for precise 
definitions of these families. The proof of the theorems is rather construc- 
tive and leads to efficient algorithms that are being implemented. 

There are several physical interpretations of FK models. The original mo- 
tivation I1FK39II was dislocations in solids. In the interpretation of IIALD83II 
Ui are the positions of a deposited material over a substratum. The inter- 
action of the atoms with the substratum is modeled by the term V. The 
periodicity of V considered in [IALD83II corresponds to a periodic substra- 
tum (e.g. a crystal) and the quasi-periodic models considered here could 
appear in cleaved faces of crystals or in quasi-crystals. 

In the interpretation of deposition, the existence of quasiperiodic solu- 
tions implies the existence of a continuum of equilibria, so that the system 
can slide. In contrast, if the KAM tori are not present, the system is pinned. 
There have been numerical explorations of these issues in [|vE99[ [vEFRJ99[ 
IvEFJO 1 1 |vEF02l IRJ971 IRJ99II . In particular, the above references pay spe- 
cial attention to the boundary of the set of parameters for which there is an 
analytic solution (breakdown of analyticity, Aubry transition), which corre- 
sponds physically to the boundary between sliding and pinning. 

Note that the physical argument to obtain sliding only requires a contin- 
uous family of solutions, whereas it has been found that often the boundary 
is described by the breakdown of an analytic family. In Section [6^ we will 
show that all the families which are smooth enough are indeed analytic. 
From the mathematical point of view this leaves open the possibility that 
there are regimes where the solutions have only a rather low regularity. In 
the case of twist mappings, using renormalization group, this regime has 
found to exist but being a codimension one surface [|Koc08|| . 

In the mathematical literature, quasi-periodic Frenkel-Kontorova mod- 
els have been considered in [|GGP06[|AP10L which use mainly topological 
methods to study the existence of orbits with rotation number. In the peri- 
odic case, the critical points of the energy, i.e. the configurations solving Q 
can be identified with orbits of a dynamical system on the annulus T x R. 
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In the quasi-periodic case, no such identification is easy. [IGGP06 II shows 
that the quasi-periodic case can be considered as a dynamical system on a 
Cantor set ( Delone set). 

To the best of our knowledge, there is no systematic Aubry-Mather theory 
analogue to that of the periodic 1-D Frenkel-Kontorova systems (existence 
of minimizing Aubry-Mather well ordered minimizers). It seems possible 
that one could get analogues of the theory of minimizing invariant measures. 
See IIBur87[|Bur88llBur90L 

1.1. Hull function. We will be interested in solutions of (|2l) given by a hull 
function 

Un = h(n ■ o)) 

where a» g R and 

h(e) = 9 + h{e) 

with 

h{e) = ^^e^"'^-^" 

equivalently, 

h{e) = h{ea) 

where : R is a function 

keZ'' 

We denote the set of h of this type by QP(a). Later on, we always use the 
notation cr = 6a for variables in T''. 
Then (O is equivalent to 

(3) hie + co) + hie -co)- IhiO) + daVia ■ hiO)) = 
where daV = ia ■ V)V. We write (|3]) in terms of h which is 

(4) hia + oja) + hicr - coa) - Ihicr) + daVicr + a ■ hicr)) = 

1 .2. External forces. We will find it convenient to add an external param- 
eter to the equilibrium equation (jU and to generalize to situations where the 
forces do not derive from a potential. So we will be looking for solutions of 

(5) &[h, A]ie) = hicr + coa) + hicr - coa) - Ihicr) + Uicr + a ■ hicr)) + A = Q 

where U : ^ M. and ti is as before, /I e R. Note that, both h and A are 
unknown. 

Adding the extra term A will allow us to study the equation ^ for ar- 
bitrary periodic U even if U is not the gradient of some periodic function. 
Later, a very simple argument developed in Section ISlLemma l22l will show 
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that, when U has a variational structure, then A = Q. Hence, in the case that 
U = daV, the equation ^ is equivalent to (O. 

Remark 1. This procedure of adding an extra parameter and showing that it 
vanishes when there are geometric properties similar to the extra parameter 
method introduced in [|Mos67ll to prove "translated curve theorems" and 
developed later in [lYoc921 |Sev99i The main advantage is that the extra 
term allows a more efficient iterative procedure. The proofs of existence of 
perturbative expansions are also considerably easier 

Our main results will be "a-posteriori" theorems (see Theorem [TT] and [T2l 
for more details) which show that if we are given a pair [/?, A] that solves 
dl]) very approximately (and provided oj, a satisfy some appropriate Dio- 
phantine condition and \h. A] satisfy some non-degeneracy condition), then 
there is a true solution close by. We note that we do not need that the system 
is "close to integrable", we only need that there is an approximate solution 
of the equilibrium equations. Of course, in the case that the system is close 
to integrable, the solutions in the integrable case are approximate solutions. 
One can, however obtain approximate solutions in other cases, for example, 
the result can validate numerically produced approximate solutions. 

The method of proof is based on an iterative procedure of Nash-Moser 
type using the quadratic convergence to overcome the small divisors. The 
quadratic convergence will be based on some geometric cancellations that 
come from the variational structure of the problem. 

We note that the method of proof leads to very efficient algorithms, which 
we present in Section 14.21 and which are implemented in [IBHdlLllll . We 
note that the method is based on the Lagrangian proof in [|Ran87[ ILMOlll 
and in [dlLOl, Section 5.3]. In [IdlLOSH . it is shown that, in the classical 
case, when V is periodic, the method of proof extends to interactions which 
are not nearest neighbors. We hope that such an extension is possible in the 
quasi-periodic case considered here. This extension to interactions which 
are not nearest neighbors is interesting because the most physical models 
involve such interactions. 

The precise formulation of Theorem [TT] and Theorem [12] requires us to 
make precise the sizes of functions, Diophantine condition, etc. which we 
take up in Section |2T[ 

Note that the solution of © are not unique. If {h{a). A] is a solution, for 
any /3 6 1., \h{a + /3a) + /3, A] is a solution. Hence, by choosing JS, we can 
always choose our solution normalized in such a way that 

1 r% . r . . „ 



(6) 




Note that the choice of /3 that accomplishes Q is unique. 
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We will indeed establish that the solution of ^ and © is unique. 

2. Function spaces and preliminary estimates 

In Section im we collect several standard definitions of spaces and present 
some preliminary results on these spaces. In Section |Z2] we present defini- 
tions of the Diophantine properties we will use in this paper. In Section [23] 
we present well known estimates for cohomology equations, which are the 
basis of the KAM procedure. 

2.1. Spaces. Given a function 

keZ'' 

it is natural to consider the function 

keZ'' 

where cr e C'/Z^. 

Clearly, if the above sums converge pointwise (in our applications they 
will converge in much stronger senses), we have 

(7) kO) = ka6) 

We will find it more convenient to define spaces and norms using the func- 
tion h. 

Following HCdlLlOl we will find it convenient to use at the same time 
spaces of analytic functions and Sobolev spaces. Using the abstract results 
in [jCdlLlOII this leads automatically to several interesting corollaries such 
as bootstrap of regularity of solutions and a numerically verifiable criterion 
for the breakdown of analytic solutions. See Section |6l 

Definition 2. If the series kv) = Z^Z'' k^^'^^'^^ defines an analytic function 
on Dp = { 7] \ \Im{7]i)\ < p} which extends continuously to Dp, we denote 

(8) ||^||p = su_ptol 

TjeDp 

We denote by s^p the space of such functions. As it is well known, £^p is a 
Banach space when endowed with ([8]). Clearly for p' < p, ^ c £/p>. By the 
maximum principle, \\f\\p is monotone increasing in pfor f e 

For the sake of convenience we will also introduce s^p, the Banach space 
of functions whose r derivative is in £/p, endowed with the norm \\f\\j2^^ 
I) I + \{f)\ where (/) denotes the average of f. By Cauchy bounds, if 
p' < p, then sz/^, c 
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We also define 

(9) \\h\\lr = J]\k\\i + \k\^y 

keZd 

and denote by the spaces of functions for which © is finite. It is well 
known that is a Banach space and indeed a Hilbert space. 

2.1.1. Some elementary properties. There are several well-known proper- 
ties of those spaces. We just note 

• Cauchy estimates for analytic functions: 

\\D%\p.s<C{l,d)-6-' -Mp 
\k\ < e-'""'-^ • \\h\\p. 

• Interpolation inequalities: 

- Analytic case: 

Lemma 3 (Hadamard 3 -circle theorem, see USteVOl "). 

. i .1 

mi < 



lip — ll"-llp_|5 \\"-\\p+g- 



- Sobolev case: 



Lemma 4 (See [IZeh75[[Ste70ll ). For anyO<n< j, < 6 < I, 
denote / = (1 - 6)n + 6j, we have for any h e H^: 

wki < ■ \\h\c ■ mfj. 

• Banach algebra properties: 

- Analytic case: 

\/g,hes^p: |||-^||p<||^||p-||/i||p. 

- Sobolev case (see [|Ada75 !l): Let m > ^, there exists a constant 
K depending only on m, d such that for any u, v e H"\ u -v e 
H'" and we have 

||m • Vllff- <K-\\u\\h"' ■ \\v\\h"-- 

• Composition properties: 

- Analytic case (see ndlLO99[ldlL08ll '): Leth e £/p and Q c 
be a compact set. Take l = distiC' - Q., (Id + a ■ h){Dp)). 

Lemma 5. Let U : Q Che an analytic function \\U{z)\\l'^(q.) ^ 
M. Define the operator <I) acting on analytic functions by {<S>\h\){cr) 
U{o- + a- h{(r)). 

We present sufficient conditions that ensure that the operator is 
well defined and differentiable. 
* If\\h* - tiWp < L, then O^. G 
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* If\\h* -h\\p < ^, then 
(DO[^]A)(cr) = djjicr + ah{(r))K{cr), 

iio[^*] - oM - Dmw - h)\i < c ■ w - hi- 

- Sobolev case (see [ |Tay971ICdlLT0]l ): 

The following result is a consequence of the Gagliardo-Nirenberg 
inequalities. 

Lemma 6. Let f e C" and assume /(O) = 0. Then, for u e 

H'"nL°° 

\\m\\H'" < K2(\\uh^)(l + \\u\\h...\ 

where KjiA) = sup|^|<^ |£>^/(jc)|. 

In the case that m> ^, if f 6 C™"*^^, we have that 

(10) \\fo{u + v)-fou-Dfou-v\\H^ < Qm(llw|lL~)(l+l|M||//'«)ll/llc-^l|v||^.. 
The reason for ([TU] ) is that we have pointwise 
/ o (m + v){x) - fo u{x) -Dfo u(x) ■ vix) = 



We obtain the desired result using that, by Gagliardo-Nirenberg 
D^f o (u + tsv) 6 H'" and that the H'" norm is a Banach algebra 
under multiplication. Therefore, we can estimate the H'" norm of 
the integrand independently of s, t. 

Note that Lemma [6] also gives a formula for the derivative of O 
as the composition with the derivative of U. It follows that if U e 
C'"+2+^, then, O e C^^K 

Remark 7. In nZeh75ll it is shown that the interpolation inequalities are a 
consequence of the existence of smoothing operators, which play an impor- 
tant role in the abstract formulation of KAM theorem. 

We also note that Lemma [5] is rather elementary. It suffices to show that 
there are Taylor estimates with uniform constants. (As shown in [|dlL099ll , 
it is important that in the domain of U, among any two points xuXi, it 
is possible to choose a path y such that £{y), its length satisfies £(y) < 

C\Xi - X2\.) 

Actually, one can prove something stronger, namely that the operator O 
is an analytic operator when U is analytic, but we will not need as much. 

The following lemma is standard in the theory of quasi-periodic func- 
tions. 
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Lemma 8. Let a be irrational. Assume that YikeZ'' \hk\ < oo. (So that, by 
Weierstrass M-test, the series h{6) = YikeX'' hj.e^'"'''"^ converges uniformly 
over the real line.) Then, 



hn 



1 

im— h(e)de. 



lim 

T 



Proof. Given e > 0, we can find A'^ such that 

El* 



\I<\>N 



6 

< -. 

~ 3 



Hence, for all T > 



5^ J^E 

\k\>N 



de 



e 

< -. 

~ 3 



Furthermore, 



(11) 



1 



hke^''"'"^de = ho+ Yj 

0<\k\<N 



^2mk-aT ^~2mk-aT 

Inik ■ a - IT 



We see that, for all T sufficiently large, we can assume that the term in the 
sum in (fTT|) (it is a finite sum) is smaller than |. This ends the proof of 
Lemma [8l □ 

2.2. Diophantine properties. Given a e R'^ such that 

(12) \a-k\>^\k\-\ VA:eZ^-{0} 

where \k\ = + 1^2! + ••• + \kdl we are interested in the numbers a; 6 R such 
that 



(13) 



\a>a ■k-n\> v\k\ ^, Vfc 6 



{0}, neZ. 



That is, we are interested in the a» e R such that coa is a Diophantine vector 
in the standard sense. Here //, v, v, r are positive numbers. 

We will denote the set of to satisfying (fT3] ) as ^(v, r; a). We also denote 
^(r;a) = U,>o^(v, r; or). 

The abundance of Diophantine numbers in subspaces is a subject of great 
current interest in number theory rKleOTllKleOS I. Nevertheless in our case, 
it suffices the following elementary result Lemma |9] (which is weaker with 
respect to the exponents obtained than the results of flKleOlllKleOSH ). 

(t; a) is of full 



Lemma 9. If a e W' satisfies (fT2l) and t > d + v, then 
Lebesgue measure. 
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Proof. Let A > 0. Consider the sets 

n v\k\~'^ 

^u,n = {oj\ \(oa-k-n\< v\kr} = {aj\ \co -| < -} 

a • k \a • k\ 

and 

^k,nA = ^k,nn[A,\mA]. 

Clearly, 

[A, IMA] \ S!{v,T;a) = yJk,n^k,nA 

So Lemma [9] will be proved when we show 
We clearly have that ^yt,n is an interval of length 

, 2v\kr 



We also observe that ^k,nA = unless A + ^ < ^ < I.OIA + ^ or 



v\kr 



\k-a\ 

\k-cY\ k-a i-v-'j-^i 1 i^^^i 



< 7^ < I.OIA - Also clearly, 



\k-a\ k-a ' \k-a\ 

UM ^k,n,A ^ 0} < 0.02A -Ik -a] + 2 

Hence 

\Uk,n'^k,n,A\< y ^^<y ^^■(0.02A-|;t-al + 2) 
■4—' \k ■ a\ ^ \k ■ a\ 



k,n V 1^' 



<0.02A • 2v • \k\-^ + — 2 ''^'"^^"''^ <yC 



k ^ k 



□ 



2.3. Cohomology equations. To prove our results, as it is standard in KAM 
theory, we have to study equations of the form 

(14) ~^{e^ui)-~m = m) 

where r] e QP(a), a» 6 R are given and ^ is the unknown. 

We note that if we use the function ^(aO) = ^(6) we see that (fT?l) is 
equivalent to 

(15) 0((r + coa) - ^{a) = f}(cr). 

The operator solving these equations is unbounded, but it satisfies some 
"tame" estimates from one space to another that can be overcome by a 
quadratically convergent algorithm. 

Clearly, a necessary solution for the existence of solutions of (fT5l ) is 



fi{cr)dcr = 0. 
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The equation (fT5l) has been considered in KAM theory. The optimal 
results for analytic functions were proved in nRus751 and we just reproduce 
the results adapted to our notation. The results for Sobolev regularity are 
very easy. 

Lemma 10. Let r) e s^p (resp. H'^, r > t) be such that 




f}(cr)dcr = 0. 



Then, there exists a unique solution of (fTSl ) which satisfies 

(16) f ^{(T)da = 0. 

Jj'' 

This solution satisfies for any p' < p 

(17) \\Up'<c{d,T)-v-'-{fi-p'r\\mp 

( resp. 

(18) ll0ll//-<C-v-i-||?7ll//., T<s<r). 

Furthermore, any distribution solution of (fT4l) differs from the solution 
claimed before by a constant. 

The method of proof is standard. We take the Fourier coefficients and see 
that (fT4l) is equivalent to: 

0^(e2..1... - 1) = flk. 

When kitO^we can use that {e^""'-'"^ _ i |-i < c • dist(/t • aco, Z)-\ From this, 
the result for Sobolev space follows rather straightforwardly. The result 
for analytic functions - with the optimal exponent quoted above - requires 
somewhat more elaborate arguments that use not only the upper bounds in 
the denominators, but also that they are not saturated very often. We refer 
to [|Rus75ll for the original proof and to [IdlLOll for a more pedagogical 
exposition. 

3. Statement OF THE MAIN RESULTS 
3.1. Statement of the analytic result. 

Theorem 11. Leth = Id + h where h{e) = h{a ■ 9) = Y.keZ'i k ■ e^"'-^'"^ with 
ho = 0, h e and a 6 R.'' such that a ■ i Q, i e - {0}. Denote 
/ = 1 + dji and T-a,a(cr) = cr - oja. We assume 

(HI) Diophantine properties (fT3l) ." \coa-k-n\ > v\k\''^, Vfc e n e 

Z. 
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(H2) Non-degeneracy condition: ||/(cr)||p < A^+, ||(/(cr)) '||p < A'' and 
\{jY-^ — )| > c for some positive constant c where (/) denotes the 
average of the periodic function f. 

(H3) \\S[h, A]\\p < e where S is defined in ©. 

(H4) Composition condition: Take l = distiC^ -Q., {Id + a-h){Dp)) where 
Q. is the domain ofU. We assume that \\h\\p + p < ^i. 

Assume furthermore that e < e*(N~ ,N'^ ,d,T,c,L,\\U\\c2(a)) ■ • p^'^*^ where 
6* > is a function and A e ( we will make explicit e* and A along the 
proof). 

Then, there exists a periodic function h* and A* eM. such that 

G[h*,A*] = 0. 

Moreover, 

\\h-h%<c-v-Y''-^-\m,mp, 

\A-r\<C-\\6[h,A]\\p. 

The solution [h*. A*] is the only solution ofS[h*,A*] = with zero average 
for h* in a ball centered at h in £/3p, i.e. [h*. A*] is the unique solution in the 
set 

\ge £^3, I (g) = 0, 11^ - h\\3, < —z ^ \ 

\ - - 2C{N-,N\d,T,c,C)] 

where C will be made explicit along the proof. 

3 .2. Statement of the Sobolev result. 

Theorem 12. Letm>\ + 2T and U 6 C^'^^'^^^^ Let h = Id + h where 
h{e) = h(a ■ 6) = Y.kez" h ■ e^"'-^'"^ with h e H'"^^^\ ho = Ofor any 6 e R 

I J ^ ^ 

and a e R such that a ■ j ^ 0, j e Z - {0}. Denote 1=1+ dji and 
T~a)a(.o-) = cr - coa. We assume 

(HI) Diophantine properties (fT3]) ." \coa-k-n\ > v\k\'''^, VA: e 7/-{0}, n e 
Z. 

(H2) Non-degeneracy condition: \\l{cr)\\H"> < N'*', \\il(.o-))~^\\H'" < N~ and 
\{tt^ — )| > cfor some positive constant c. 

(H3) \\S[kA]\\H^<e. 

Assume furthermore that e < e*(N~,N'^, d, r, c, ||i7|| )-v ^ where e* > 
is a function which we will make explicit along the proof. Then, there exists 
a periodic function h* e H and /i* e R such that 

6[h*,A*] =0. 
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Moreover, 

\\h-h*\\H.n-^r<C-v~~^{N^f -6, 
\A-A*\<C-e. 

The solution \h*,A*\ is the only solution of&\h*,A\ = with zero average 
for h* in a neighborhood of[h. A] in i.e. [h* , A*] is the unique solution 

in the set 

\g 6 H'"^'' I {g) = 0, \\g - h\\H,...^r < — 1 

[ 2C{N , N+, d, T, C, ||t/||c»'+34r+l) ■ Cm-4r,m+4r 

where C will be made explicit along the proof. 

Remark 13. The theorems [TTIand [I2]have the "a-posteriori" format of 
numerical analysis. 

Given a function which solves very approximately the invariance equa- 
tions, then there is a true solution nearby. The needed approximation is 
quantified in terms of the non-degeneracy condition N'^,N',c, (which in 
numerical analysis are often called condition numbers ). We note that the 
condition numbers can be computed in the approximate solution. We em- 
phasize that this formulation does not require that the system is close to 
integrable. 

Of course, the non-degeneracy conditions depend on the function h (and 
on the parameter of the domain p in the analytic case. If this can cause 
confusion - e.g. when we are performing an iterative step - we will use 
N%kp). 

4. Proof of main theorems. Theorem [TTIand Theorem [T2l 

As indicated in the introduction, the proofs of Theorems [TT] and [T^l are 
based on an iterative step that given an approximate solution of (fT9l ) will 
produce a better approximation. 

The crux of the proof is to show that, if started with a sufficiently approx- 
imate solution, the procedure converges. 

As it is very well known in KAM theory, there are arguments that es- 
tablish the convergence, provided that we show that the iterative procedure 
satisfies tame quadratic estimates. That is, that the norm of the new error is 
bounded by the square of another norm of the original error (in a smoother 
space) times a factor that depends on the "loss of regularity". There are 
several abstract theorems in this direction, one which is quite well adapted 
to the problem at hand and which we will use appears in HCdlLlOH . 

In Section l4n we will give some motivations for the iterative procedure. 
(It is a Newton method with a small modification that does not affect the 
quadratic convergence.) 
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In Section 14.21 we will formulate the iterative procedure as a succes- 
sion of elementary sub-steps. We note that these elementary sub-steps can 
be implemented by very efficient algorithms. If the functions h are dis- 
cretized using points and appropriate algorithms are used for the math- 
ematical sub-steps, the iterative procedure requires only 0(N) storage and 
0(N log N) operations. 

In Section B3l we present estimates for the iterative step. We first present 
estimates on how much it changes the function. Then, we use the Taylor 
estimates to show that the error of the improved function is tame quadratic 
in the sense of Nash-Moser theory. 

In Section 14. 4[ we will review the convergence of the procedure, which, 
as we have indicated before is rather standard, indeed in [ICdlLlOII . there is 
an abstract theorem designed to cover exactly the problems considered here. 
Nevertheless, for the sake of completeness in the analytic case, we will also 
present a very short direct proof of the convergence argument. Since this 
direct proof is so direct it also leads to good numerical values. 

In Section |431 we present several considerations that allow us to discuss 
uniqueness. 

Finally, in Section [6l we show that we can obtain several consequences 
combining the Sobolev and analytic versions of the a-posteriori theorem. 
Namely, we show that for analytic mappings all sufficiently smooth solu- 
tions are analytic, that Lindstedt series converge, and we present a numeri- 
cally efficient criterion for the breakdown of analyticity. As it was shown in 
[IALD83II , the breakdown of analyticity is associated to the onset of trans- 
port properties, so there is some interest in its computation. 

Of course, readers interested only in rigorous proofs can safely skip Sec- 
tion 14.11 and the algorithmic considerations in Section 14.21 Readers only 
interested in algorithms can skip 14.31 

4.1. Motivation for the iterative step. We start from an approximate so- 
lution of ©, 

(19) h{(T + a>a) + h(cr - too) - 2h{cr) + U(cr + a ■ h{o-)) + A = e 

where e is to be thought of as "small". 

Our goal is to devise a procedure that gives a much more approximate 
solution. For the moment, we will not make precise the sense in which 
quantities are small. This will be taken up in Section 14.31 

Given an approximate solution as in ([T3, the Newton method would con- 
sist in finding a solution of 

(20) A(cr + coa) + A(cr - cua) - 2A(cr) + U{cr + a ■ h{a)) ■ A(cr) + 6 = -e. 
Then [ti -I- A, A + 6] will be a better approximate solution. 
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The equation is hard to study because of the term daU((T + a- h{cr)) ■ 
A(cr) which is not constant coefficients. 

The key observation is that, if we are given (fT9] ) we are also given the 
following equation which is just obtained by taking the derivative of (fT9l ) 
with respect to 6 (we recall that a = 6a): 
(21) 

dah(cr+a)a)+dah(cr-a)a)-2dah(cr)+daU(cr+a-h{cr))-{l+dah{(T)) = e (6). 

Denoting Z(cr) = 1 + dji{a), we rewrite (|2TI ): 
(22) /V + oja) + Z(cr - oja) - 2l{a) + djj{a + a- h(cr)) ■ tier) = e'iO). 



If we substitute (|22|) into (|20l) . we obtain that the Newton procedure (|20l ) 
is equivalent to 

(23) A(cr+a;a)+A(o--aja)+^^ ^ ^ ^■A((r) = -e-6. 

1{(t) 

The equation (l23l) is still hard to solve, but, as we will now see, we can 
study the equation obtained by omitting the term e'{6) ■ A(cr) in (1231) . 

Because both e'(6) and A((t) are small, we can hope and we will show 
that omitting the term e' ■ A((t) does not affect the quadratic character of the 
procedure. 

The key observation is that the quasi-Newton equation for A, 6: 

(24) _ 

l((r)-A(cr+oja)+l{o-)-A(cr-a)a)-A(cr)[t((r+ioa)+t((r-ioa)] = {-e-6)-l{a) 
is equivalent to the system 



(25) 



A^ 






J. 




J. 





w 



l-loT. 



(26) W oT^c^-W = l-{e + 6). 

In fact, from (|25] ), we get 

w = Aor_^„./-A-/or_^„ 

and 

o = A ■ z o - A o ■ / 

By (|26l) . we can easily get the equivalence. 

The key point is that the equations (|25l) and (|26l) are of the form (fT4l) . and 
can be studied using the theory in Section 12.31 

We write W = + W where H^" is a function with zero average and 

V7 is a number. That is, we decompose W into its average and the zero 
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average part. Both are unknowns. Now we describe the procedure to solve 
the system ([25]), (l26l) . 

We first choose 5 to be the unique value that makes the average of the 
right-hand-side of (1261 ) zero. Then, we can apply Lemma [TO] to find 
solving (|26|) . We note that there is only one choice of 6 and then, is 
determined uniquely, by the condition that it solves (1261 ) and that it has zero 
average. The only solutions of (|26l ) differ from it by a constant. 

Then we observe that W = - 1 "° | is the only possible value of the 



average of solutions of (|26l ) that makes the right-hand- side of (1251) with 
zero average. Them, we can apply again Lemma[lO]to find j solving (l25l) . 
This solution is unique up to the addition of constant. Once we have j, we 
obtain A is obtained just multiplying by /. Note that the A is thus determined 
uniquely up to the addition of a constant multiple of /. In particular, A is 
unique when we impose the normalization that it has zero average. 

4.2. Formulation of the iterative step. 

Algorithm 14. Given h : T'' ^ R, A £ R with h{o-) = Y.kei" he^""'"' and 
h(6) = h{aQ) for e R and any irrational vector a e R , d > 2, we will 
calculate: 

1) Let ^ = h{cr + aco) + h{a - aoS) - lh{a). In Fourier components 
S£k = 2(cos (jja ■ k - l)hk. 

2) We calculate U(cr + a ■ h{cr)). 

3) So we can calculate e = ^ + U + A. 

4) Calculate 1=1+ dji. In Fourier components 4 = 5/t,o + ■ a-hk 
where 6o^k is the Kronecker delta. 

5) Let f = !■ e. 

6) Choose 6 = -(/). 

7) Denote b = l-(e + 6). 

8) Solve the cohomology equation (|26l) for W'^ with zero average. That 

is, 2(cos (jff-^-l)' 

— <^^> 

9) Take W = -J±l^^. 

10) Calculate W = W° + W. 

11) Solve the cohomology equation (1251). Find satisfying fi-fioT-^ja = 

— . That is, Bh = — , .x where a = — . 

12) We obtain A = ■ I where is chosen to be -0 ■ I) I {I) so that 
(A) = 0. 
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Remark 15. It is important to note that the procedure also shows that the 
solution of (|25l ) and (|26l ) is unique up to the addition of a constant multiple 
of /to A. 

Hence, if we the choose the solutions of (|25] ). (|26l ) which satisfy the 
normalization (A) = 0, the solutions are unique. 

Of course, since the system (|25T ). (|26l) is equivalent to (|24l) . the same 
considerations apply to (l24l) . 

Remark 16. Note that if we consider functions discretized by their values 
at points and by A^ Fourier coefficients, steps 1), 4), 8), 11) are fast in 
Fourier coefficients (they require 0(N) operations) while the other steps are 
fast (they require 0{N) operations) in the representation of the function by 
its values at points. Of course, once we know the representation in space or 
in Fourier coefficients we can use the Fast Fourier transform which requires 
0(N log N) operations to compute the other. 

Note also that if we discretize the function as above, the iterative step 
only requires to store several functions, and therefore we only need to store 
0(N) numbers. 

We, thus obtain a quadratical convergent algorithm, with 0{N) storage 
requirements and 0{Nlog{N)) operations. In contrast, a Newton method 
would require 0{N^) storage to store a matrix and 0{N^) operations to 
solve the linear equations (there are faster algorithms [|Knu811 to solve lin- 
ear equations but they do not seem to be practical). In practice the present 
algorithm with A^ = 10^ can run comfortably on a modest desktop machine. 

Remark 17. It is important to note that [A, 6], the outcome of the algorithm 
depends linearly on e = S[h, A]. 
Hence, we will write 

(27) [A,6] = r][h,A]e 

The operator rj is called an "approximate right inverse " in Nash-Moser 
theory. See, for example [|Zeh75l . 

Notice that the estimates for the improved solution can be written in a 
symbolic way as estimating &\\h. A] + r]\h, A]&\h, A] ], which, using Taylor 
expansion (up to quadratic errors) becomes 

(28) &{K ^] + D&{h, A]ri{h, A]&{K ^] 

In the Newton method, we would choose rj in such a way that (|28] ) vanishes. 
As pointed out in [|Mos66irZeh75ll , it suffices that the norm of can be 
bounded by the square of another norm of A]. 

4.3. Estimates on the quasi-Newton step. In this section we show that 
the Quasi-Newton method specified in Algorithm [14] produces more ap- 
proximate solutions. We will present two versions of the estimates, one in 
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analytic spaces and another one in Sobolev spaces. The goal is to obtain 
that the new error is quadratic in the original error even if in a weaker norm. 
We note that the analytic estimates presented are a bit more delicate and 
involve a condition, (l34l) . 

4.3.1. Some useful identities. We start by remarking an elementary identity 
that will be used for both the analytic and the Sobolev estimates: 

(29) I ■ (Di6[h, A]A) - A • (DiB[h, A]l) = -I ■ i6[h. A] + 5). 

where Di denote the derivative with respect to the first variable. 

We also have the following identity obtained just adding and subtracting 
terms in the definition of B[h + A,A + 6] and grouping them. 

(30) 
B[h + A,A + 6] 

=6[h. A] + A(o- + a>a) + A(cr - coa) - 2A((t) + 6 

+ U{cr + a-{h + A)(cr)) -U{cr + a- h{o-)) 

Pfu ji^/ fft , kcr + coa) + tjcr - coa) - lljcr) ^ 
=6[h. A] + (-6[h, A]) + A(cr) 

Kcr) 

+ Uicr + a-ih + A)(cr)) -U(cr + a- h{o-)) 

e'(e)-daU(o- + a- kcr)) -tier) ^ ^ . _ _ 

=— ^— -^^ y y ^ . ^^^^ + [/((^ + Q, . (/j + A)((r)) -U{cr + a- h(cr)) 

l{a) 

=e' ■ 4^ + C7(cr + a • (^ + A)((r)) -U{(T + a- Her)) - daU((T + a ■ h) ■ A(cr) 
1((t) 

Kcr) 

where R = U(cr + a-(h + A)(cr)) - U((t + a ■ h{cr)) - djj{a + a-h)- A(o-). 
Clearly, R is the remainder of the Taylor estimate in the composition studied 
in Lemma [5] and Lemma [6l 

4.3.2. Estimates for the iterative step in analytic spaces. We now observe 
that for any p' < p, by (fTTl) . we obtain using (l26l) 

\\wX<C{d,T)-v-' -{p-p')-' -N^ -WeWp. 

Since the average of W is obtained in step 9) of Algorithm (fT?l) . we have 
the estimate for W: 

\W\<c- WW^Wp, ■ {N-f < c ■ {N-f ■ ■ C{d,T) ■ ■ \p-p'r ■ \\e\\p. 
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Therefore, we obtain the estimates for W: 

\\W\\,><M-v-'-(p-p'r -Well, 

where M = (c-(N')^ + l)-N'^C(d,T). The important point is that the constant 
is uniform provided h stays in a neighborhood in || • ||p norm. 
Again for p" < p', by (fTT]) . we have 

\\~A\\p" < Cid,Tyv-'-(p'-p"r-iN-f-\\W\\,, < M'-(p-p'r-(p'-p"r-\\e\\,. 
So, we have 

mip" <N^-M- v-^ip-pT ■ ip' -P'T ■ \\e\\p 

and 

(31) IIAIIp. < A^^ • N-M ■ v-^(p-pT ■ ip'-p'T ■ Ikllp. 
Similarly, using Cauchy estimates, we obtain for p'" < p" 

(32) m\pn < ■ N-M ■ v-^(p" -p"r'(p-p'r ■ (p' -pT ■ mi. 

If we take p - p' = p' - p" in (|3T1) and p - p" = p' - p" = p - p" (and 
redefine p'" ) in (l32l) we obtain: 

^^^^ IIAIIp. <M'-y-V-p"r''-|Hlp 

m\p"<M'-v-^{p-p"r^'-'-\\e\\p. 

If we have that ||A||p" < t/2, which, by (1331) is implied by 

(34) M' -v-Hp-p")-^' ■\\e\\p<il2, 

we can define U((t + ah + A(cr)) and indeed apply Taylor's estimate we 
obtain 

(35) < sup C7 • IIAIlJ, < Mv-'(p" -pr^'WeWj. 

The first term in the right-hand- side of (|30^ is estimated using the Cauchy 
estimates and the previous estimates on ||A||. 

lie' • 4|l ,„ < (p-p"r'\\e\\pN-N* ■ M' ■ y-^ -p")"'' • Ikllp 



= M(p-p")"'-iy-2|H|2 



<M(p-p"r^'-'v-'%f 



lip 

||2 

Hp- 

The last estimate is done with the purpose of simplifying the expressions, 
but it is obviously wasteful. Note that r > 1 and that the estimates above 
are delicate only when p - p", y are small. 
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Finally, putting together the estimates for the two terms in the right-hand- 
side of we have: 

(36) \m + A, i + 6%,, < C ■ V'ip - p'^wm, 
Therefore, we have proved the following 

Lemma 18. In the hypothesis of Theorem [77] 

Assume that (l34l) holds. Then, the improved function obtained applying 
Algorithm (fT4l) . satisfies (|36l) . 

As it is well-known in KAM theory, the above estimates imply that the 
iterative procedure can be repeated indefinitely and the resulting sequence 
converges to a function satisfying the claims of Theorem [TT] Indeed, the 
paper HCdlLl O] contains an abstract theorem that immediately applies to 
this situation. We will discuss this in more detail in Section 

4.3.3. Sobolev estimates for the iterative step. Let s > j. According to the 
algorithm 8), we have 

\\b\y = J](l + \k\y\k-(e, + 6t,o)\' 

keZd 

= Yj (1 + l-tPr • 14 • 

<\\f\y<K-\\i\y-\\e\y. 

By we obtain 

\\W%s-r < C • y-^ ■ IMhs <Cv-^ KN^ ■ llelb... 

We get 

\W\ < C ■ iN~)W\y-r 

and 

\\W\\H.^-r <C v ^ -K-N^ -WellH^. 

So we have 

||A||^.-2. < C ■ v-'K^ ■ (N-f ■ \\W\y-r <C-v-^-N^- (N-f ■ \\e\y. 
Hence, 

\\A\\H^-2r<C-v-^-(N^f-iN-f-\\e\y. 
We recall that the approximate inverse of the derivative ri[h. A] is just the 
result of applying applying the algorithm in Section 14.21 i.e. [A, 6] = 
ri[h, A]S[h, A]. We have proved the following lemma: 

Lemma 19. Let s > j + 2t. Then we have 

\\i][h,A]6[h,A]\\H..-2r < C • • (N^f ■ (N-f ■ \\e\\H^. 
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We will also need estimates on {DiS[h, A]rj[h, A] - Id)(S[h, A] + 6). 
Lemma 20. 

\\iDi6[h, A]7][h, A] - Id)i6[h, A] + 6)\\H^-2r 
<C ■ • (N^f ■ (Ny ■ mk ^]||^.-2.-,||fi[^, ^]||^.. 

Proof. By the definition of rj[h,A], we know A = -rjlh, A](S[h, A] + 6). 
Hence, 

iDi6[h, A]r][h, A] - Id)iB[h, A] + 6) 

. - . 'A-Di6[h,A] 
= - DiS[h, A]A - &[h. A]- 6= . . 

So we have that 

\\iDiB[h, A]r][h, A] - Id)iB[h, A] + 6)\\H^-2r 
<C ■ v'^ ■ (N^ ■ (Ny ■ \\6[h,A]\\H.^-2r-mkA]\y. 

□ 

4.4. Convergence of the procedure. The existence of solutions both in 
the analytic case and in the Sobolev case is deduced from the estimates in 
Section |43] followed by Nash-Moser estimates. 

Indeed in HCdlLlOL one can find an abstract Nash-Moser implicit func- 
tion theorem which is tailored to the theorems [TTIand [T2l 

In this section, we reproduce the theorem from r CdlLlOl and explain why 
it is applicable. We note that the the theorem has several corollaries which 
are of physical interest and we present them in Section |6] 

For the sake of completeness, in Section |4. 4. 2[ we present a direct proof 
of the convergence in the analytic case. 

4.4. 1 . An abstract implicit function theorem. In HCdlLlOl Appendix A] one 
can find a proof of the following result. Theorem |2T1 This is an abstract 
theorem that applies to operators in scales of Banach spaces, which have 
smoothing operators. 

In [Cdl LlOU one can also find a verification that the Sobolev spaces and 
analytic spaces considered indeed have smoothing operators (one can take 
S tTik^k^^'"'"^ = e~'^''^hi,e^^'''"' . The regularity properties of the operator en- 
tering in the assumptions of Theorem IITI follow immediately for the com- 
position properties presented in Section 12. 1.1 1 and, specially Lemma[6l 

Theorem 21. Let m > 2t and X'' for m < r < m + 34t be a scale of Banach 
spaces with smoothing operators as shown in HCdlLlOH . Let Sr be the unit 
ball in , Sr = h + S,. the unit ball translated by h £ X'' and S{X'', X''^^'^) 
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is the space of bounded linear operators from to Consider a map 

T -.Br^ X''^^ and rj : Sr ^ S{X\X'-^') satisfying the following: 

(i) T{Sr n X') c X'^^^for m<r<m + 3 At. 

(ii) "Fls^nxr : Sr C) X'' ^ X'"^'^ has two continuous Frechet derivatives, 
both bounded by some constant M, for m < r < m + 34t. 

(iii) ||77[A]:r[A]^-2. < C • iif [A]|U- , A 6 SrJor r = m - 2t, m + 32t. 

(iv) WiDTiAMA] - Id)r[A]h--2r < CJ|T[A]||^,., A 6 !Br,forr = m. 

(V) ||r[A]|U,„.32. < C • (1 + ||A||^,„.340, A 6 S™. 

Then if \\'F[h]\\x'"-2r is sufficiently small, there exists h* e X'" such that 
r[h*] = 0. Moreover, \\h - < C • ||rMIUm-2r. 



We recall that the method of proof of Theorem 121] following [|Sch60ll 
is to modify the quasi-Newton step adding a smoothing step. That is, one 
constructs a sequence [/?„+i, /l„+i] = [h„,-^n] + S t„ri[h„, An]&[hn, ^n]- The 
choices of have to be carefully chosen so that the quadratic convergence 
(in some norm) is maintained. The main difference between Theorem |2T| 
and the result in [lSch60ll is that Theorem IITI includes the fact that 77 is an 
approximate inverse and not an inverse. 

The estimates showing that rj is indeed an approximate inverse are the 
estimates obtained in Section |431 



4.4.2. A direct proof of the convergence in the analytic case. Since the es- 
timates in the analytic case are so easy, we present a direct proof. As we 
will see, the estimates are rather easy to verify. The main difficulty is the 
order of the choices. 

We start with an approximate solution [^0, -^o] with h 6 .s?/^^. 

Since we will have to change the function through an iterative procedure, 
we note that the condition numbers N'^,N~,c depend on the functions we 
are considering, nevertheless, they are uniform in a neighborhood. 

We start by choosing a number y > such that that in neighborhood of 
size y in £^[, we have that A^* < IN-Qiq), c < IcQiq). 
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The following algebraic identities will be useful in estimating the change 
of non-degeneracy conditions in the iterative step. 

NHh-p,,) = \\\+dM\p, 

<N^(ho;po)\\da(h-ho\\po 

< N-{ho;po) + \\da(h - ho)\UN-(ho;po)N-(h;p) 



1-1° T-cja Iq • lo ° T-cja I 
/q/o O T.^Jl O T-^a I 

< [a^-(^;po)A^-(^o;po)]'(II/IU + ll^olU)!!^^- ^^oll, 



Hence we can find a number y > depending only on the non-degeneracy 
conditions N-,c so that all the functions in a ball of radius y in £/p^ cen- 
tered at ^0, have non-degeneracy constants not larger than twice the non- 
degeneracy assumptions of fiQ. 

More generally, we have, by the same argument that if \\h - hoW^^t < y, 

then, N^ikp) < 2N-*{ho,p). 

The key estimates are, as follows to show that, with some convenient 
choices of radii, which we do at the outset, the iterative process can be 
applied indefinitely and indeed it converges. We will use (|37]) to show that 
the non-degeneracy constants do not deteriorate much. 

We denote by 



(38) Pn = Pn-l - j2-" = Pod 



and provided that we can apply the iterative step (that is, provided that (l34l) 
applies with the choices of p„ in (|38l) . we define for n > 1, [h„,A„] = 

[hn-l,A„-i] + ri[hn-i, An-l]B[h„-i, An-l]. 

We denote by M the constant in Lemma [TS] corresponding to twice the 
degeneracy assumptions corresponding to the original function. 



(37) 



\cih) - ciho)\ = 
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If (|34l) applies n times, for typographical simplicity, we denote 6,- = 
\\&[hi, Mlpi we see that 



(39) 



2^ 

2 



< (^y-2p-4ry+2+-+2"2(n-l)4r+2(n-2)4T+- 



2"- Mr ^2" 
^0 



We see that if {Mv~^p~'^^)^2^''eQ < 1, the right-hand- side of (|39] ) decreases 
faster than any exponential. Indeed the factor can be made as small as 
desired by assuming that 6o is small enough. 

If we apply n-times the inductive step, we see that the distance from the 
range of /i„ to the complement of the domain of definition of U is at least 

n n n 

(=0 i=0 i=0 



i=0 

Note that if eo is small enough, this is bounded from below by i| indepen- 
dent of n. 

According to Lemma [TSl the only thing we have to verify is (|34] ). which 
with the choices of radii that we have made amounts to: 

We note that this condition is satisfied independently of n if n is large 
enough. 

Using (|37]) . we have: 

N\k,Pn) < N\K.uPn) + \m„\\p„ < N\K.,,Pn-l) + Mv'^ pl'^'2''^"~'^ Sn-l 

m-1 



i=0 

and similarly for N~,c. Therefore, under smallness conditions on 6o, we 
get that the non-degeneracy conditions do not change by a factor 2 from the 
original one, so that the induction hypothesis are satisfied. 

In summary, under just three smallness conditions in 6o, which can be 
assessed just looking at the non-degeneracy conditions, we conclude that 
the iterative step can be carried out infinitely often and that the assumptions 
on the non-degeneracy constants make in the estimates for the step remain 
valid. 
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We also note that since p„ > po/2, we have 

N 

\\hN - ho\\po/2 + U« - /lol < ^ \\hn - hn-l\\po/2 + Un " ^n-ll 

n=l 

N 
n=\ 

which establishes the quantitative claims made for the result. 

4.5. Uniqueness of the solution. In this section, we establish the unique- 
ness claims for the Theorems [TT] [l2l We note that the proof is very el- 
ementary and only uses the theory of linearized solutions as well as the 
interpolation inequalities in Section 12.1.11 

4.5.1. Uniqueness for the analytic case. If \\k* — ^**||£,|/l* - X'*\ is suf- 
ficiently small and £[^*,i*] = fi[A",i"] = 0, by Taylor's theorem and 
Lemma [51 we have 

(40) = fi[r*, /I"] - X*\ = Di6[h\ A*](h** - h*) + {A** - A*) + R 
where \\R\\p < C ■ \\h**-h*\\l. 

Now, denoting as before 1=1+ dJi* and recalling that 

Di&[k,A*]-l= ^&[k,r] =0 
0.9 

we can write the equation (1401 ) as: 

(41) / • {Dy&[k,A*]{k* - k)) - (k* - k) ■ {Di&[k,A*\i) = -ir. 

The proof of uniqueness is based on uniqueness of the solution of the system 
(l25l) and (l26l) . By the estimates in Section [431 we conclude that for any 
0<p"<p'<^, 

11^" - k\\p. < ciN-,N\d,T,c) ■ • (p-pT ■ ip' -p'TWWp 

Takep" = f andp' = ^p. 

In the analytic case, we obtain 

11^" -k\\p<c- • p-2^ • 11^" - k\\l 

8 4 

<c ■ • p-^' ■ \\k* - k\\p ■ \\k* - kh., 

8 8' 

where C > is a constant depending on N~ ,d,T,c,C. The last in- 
equality holds by Lemma [3l So when \\k - h\\ip small enough, we obtain 

k* = k. A** = A*. This completes the proof of uniqueness of the solution 
in Theorem [TT] for the analytic case. 
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4.5.2. Uniqueness for the Sobolev case. Instead of applying Hadamard 3- 
circle theorem for the analytic case, we use the interpolation inequality for 
Sobolev case ( Lemma U). 

Following the proof in Section IT. 5. 11 we will have 

<C V' ■ Cn,~AT,m+4T ' \\h** - /z*||//m-4. • \\h** - /z*||//m+4r. 

This completes the proof of uniqueness of the solution in Theorem [T2] for 
the Sobolev case. 

5. VANISHING LEMMA 

In this section we prove 

Lemma 22. Consider a solution of Q with the stated periodicity condition. 
If 

U = daV 

then A = 0. 

Proof. The proof is very simple. We multiply (|5]) by h'{6) and compute 
limj-^oo 27 J_7 of the terms. We note that this produces the formula 

1 r^^ 

(42) A = -\im— U{hm ■ h'me. 

In fact, we observe that 

h{e + oS) + hie -co)- 2h(e) = h(e + co) + Mo - - 2h(e) e QP{a). 
Similarly, 

h'{e) = 1 + h'{e) e QP{a). 

Hence, 

W + oj) + h(e -co)- 2h{e)] ■ h'iO) 6 QP{a) 
and we have that 
1 

lim — \h{e + a>) + h{6 -co)- 2h{e)] ■ h'(e) 

r-»oo 2T j^j 

= ^ -hk ■ 2{cos(2ncok ■ a) - I) ■ h-k2m(k ■ a) + ho{2 cos{2nco0 - a) - 1) 

keZ''-[0\ 

= 0. 

The first equality is true because of Lemma [8] and the fact that the sum is 
Cauchy formula for the k = coefficient of the integrand ( as we will see 
below). The fact that the sum is is clear because it is antisymmetric in k. 
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In fact, note that 

[h o + h o - 2h]i^ = 2(cos {Incok • or) - \)hk, 

in particular, the coefficient vanishes for k = Q, and 

\h'\ = 6Q^k + ■ a ■ hf,. 

We have using Cauchy formula for the Fourier series of the product 

[[hoT^ + h7T.^-2h]-h% 
= y 2(cos(2;ra»fc • a) - l)/z/t • [5o,jt ~ 27r/A: • a • 

keZ^ 

= ^ -^i- ■ 2(cos(27ra»fc ■ a) - 1) • h.i,2ni(k ■ a) + ho(2 cos(2nco0 ■ a) - 1). 

keZ''-{0] 

We also observe that 

J d^ViahiO)) ■ h'{e)de = V{ah{T)) - V{ah{-T)). 

So it is bounded independent of T. When we divide the integral by 2T and 
take the limit T ^ oo. We obtain 0. This ends the proof of Lemma l22l □ 



6. Several further consequences of the formalism 

As pointed out in HCdlLlOI . once one has an a-posteriori theorem with 
local uniqueness in analytic and Sobolev spaces, there are more or less auto- 
matically several consequences which could be of interest for applications 
and which we now make explicit for our case. 

6.1. Existence of perturbative expansions to all orders and their con- 
vergence. If we consider models in which the interaction has a small pa- 
rameter, i.e. the interaction term is given by eU, it is interesting to know 
whether one can write formal power series = 2„ e"h'\ = X« ^"^^'^ 
which solve ^ in the sense of power series as well as the normalization 
condition Q. Furthermore it is interesting to show that that the series con- 
verges. These power series for hull functions are very similar to the Lindst- 
edt series in mechanics. 

We will show that, when the frequencies are Diophantine, the solution to 
both questions is affirmative. Series exist to all orders and converge. 
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6.1.1. Existence of Lindstedt series to all orders. We first argue that one 
can find the solution to © in the sense of power series. 

If we substitute the power series and match like powers of e, we obtain a 
hierarchy of equations for the coefficients of the perturbation. At order 6° 
we obtain: 

(43) h^io- + coa) + h\cr - coa) - 2h\(r) + = 

which implies that /l" = 0, is a constant. Because of the normalization 
we have = 0. 
At order , we obtain 

(44) h\o- + coa) + h\o- - coa) - 2h\o-) + U{o-) + =0. 

This equation is very similar to the equations studied in Section [23l Indeed, 
in Fourier series, it is equivalent to 

hl2(cos(2na)k ■a-l) = Uk + 6o^k^\ 

We see that we can determine = -Uq. is not determined by (HD but 
the normalization <^ sets it to /z^ = 0. All the other Fourier coefficients 
can be determined and indeed if Uj, is analytic in some domain then hj^ is 
analytic in a slightly smaller domain. 
In general, at order n, we obtain: 

(45) h"{o- + coa) + h"{(r - coa) - 2h"(cr) + R„{cr) + = 

where 7?„ is a polynomial expression in h\h^, ■ ■ ■ ,h"~^ with coefficients 
which are derivatives of U. So, we can assume by induction that R„ is 
known. The equation (l45l) is of the same form as (|44] ) and the same analy- 
sis shows that we can get a unique solution for h" and, hence, recover the 
hypothesis. 

6.1.2. Convergence of the formal power series. The fact that the formal 
power series converges is a very easy consequence of the fact that there are 
analytic families h^, A^, which solve the equations. This is a general fact, 
which is a consequence of the formalism and we go over the proof rather 
quickly. See PCdlLlOllGEdlLOgH . 

We just note that the method we have used works just as well for complex 
functions. We note that for all e small enough, there is a solution. (Note 
that we can take ^ = 0, /I = as an approximate solution if e is small). So, 
it suffices to show that this solution depends differentiably on the complex 
parameter e. We follow the standard practice in analysis of first obtaining a 
guess of the derivative and then proving it indeed satisfies the definition of 
derivative. 
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For a fixed e we can guess j-J^^, j^A^ because if tliey existed, they should 
satisfy 

d /y d ^ d ^ 

—h^{a + (jja) + —h^{(T - coa) - 2—h^(cr) 
(Arx de de de 

^^^^ ^ . d. ^ . d 

+ edaU{cr + ah^)—he{cr) + U{a + ah{a)) + —Ae = 0. 

de de 

The method used in Section |4?n shows that the equation (|46l) for^^£(cr) can 
be transformed into a constant coefficient equation (note that, by assumption 
is an exact solution of ([5])). 

Now, to prove that this guess indeed is the derivative, we just note that 
+ fij^he)\\p-,j < C|/zp. Then, applying the a-posteriori format and 
the local uniqueness, we conclude WK+i^ -K- A'^^ell(p-^)/2 ^ Cfj^. 

6.2. Bootstrap of regularity. In this section we state the theorem of boot- 
strap of regularity and omit the proof. See [ICdlLlOl for more details. 

Theorem 23. Let h e H"\ Abe a solution of ^ with U analytic. 

Assume that m is large enough (depending only on the Diophantine ex- 
ponent). Then, h is analytic. 

The idea of the proof is very simple. We can take a truncation of the 
Fourier series as an approximate solution. Of course, these are analytic 
functions. If the decrease of the Fourier series is fast enough, it is possible 
to use the analytic theorem and conclude that there is an analytic solution. 
By the uniqueness in Sobolev spaces, this must be the original solution. 

In Sobolev regularity this is restated as Theorem 6.8 in [ICdlLlOll . In 
[|SZ88[ iGEdlLOSH . one can find a similar argument for classes. The ar- 
gument for C classes in the later papers is somewhat more involved since 
it obtains sharper results by using better approximations than truncating. 

6.3. A practical numerical criterion for the analyticity breakdown. The 

above considerations lead to a very practical and reliable way to compute 
thresholds of breakdown of analytic solutions. 

Observe that by now, we have efficient algorithms to compute the invari- 
ant tori, given approximate solutions. This, of course, immediately leads 
to a continuation algorithm. Since we have an a-posteriori theorem, we can 
be sure that, the approximate solutions produced numerically (which satisfy 
the invariance equation up to a few units of round-off error) correspond to 
true solutions if they satisfy the non-degeneracy conditions. 

Therefore, a practical algorithm to compute the threshold of breakdown 
is to implement the continuation method and monitor the non- degeneracy 
conditions. 
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In many occasions it happens that the only condition that fails is that 
becomes very large. In that case, one can argue that the solutions ex- 
perience a breakdown because if there were analytic tori in a neighborhood 
of parameters, the Sobolev norms would remain bounded. Similar meth- 
ods for the periodic Frenkel-Kontorova models and models with long range 
interactions have been implemented in [CdlLlOJ. 

Some implementations of the method are already in progress [IBHdlLllll 
and the results will be reported elsewhere. 
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